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ABSTRACT 

We present an analysis of the broad-band spectra of 4U 1705-44 obtained with Suzaku in 2006-2008 
and by BeppoSAXva 2000. The source exhibits two distinct states: the hard state shows emission from 
1 to 150 keV, while the soft state is mostly confined to be < 40 keV. We model soft-state continuum 
spectra with two thermal components, one of which is a multicolor accretion disk and the other is a 
single-temperature blackbody to describe the boundary layer, with additional weak Comptonization 
represented by either a simple power law or the SIMPL model by Steiner et al. The hard-state 
continuum spectra are modeled by a single-temperature blackbody for the boundary layer plus strong 
Comptonization, modeled by a cutoff power law. While we are unable to draw firm conclusions about 
the physical properties of the disk in the hard state, the accretion disk in the soft state appears to 
approximately follow L oc T^'^ . The deviation from L cx T^, as expected from a constant inner disk 
radius, might be caused by a luminosity-dependent spectral hardening factor and/or real changes of 
the inner disk radius in some part of the soft state. The boundary layer apparent emission area is 
roughly constant from the hard to the soft states, with a value of about 1/11 of the neutron star 
surface. The magnetic field on the surface of the NS in 4U 1705-44 is estimated to be less than 
about 1.9 X 10^ G, assuming that the disk is truncated by the ISCO or by the neutron star surface. 
Broad relativistic Fe lines are detected in most spectra and are modeled with the diskline model. The 
strength of the Fe lines is found to correlate well with the boundary layer emission in the soft state. 
In the hard state, the Fe lines are probably due to illumination of the accretion disk by the strong 
Comptonization emission. 

Subject headings: accretion, accretion disks — stars: neutron — X-rays: binaries — X-rays: bursts 
— X-ray: stars 



1. INTRODUCTION 

There are two main classes of luminous and weakly 
magnetized neutron stars (NSs) in low-mass X-ray bi- 
naries (LMXBs), i.e., atoll and Z sources, named after 
the patterns that they trace out in X-ray color-color 
diagrams (CDs) or hardness- intensity diagrams (HIDs) 
(jHasinger fc van der Kli jfToSOt Ivan der Klisl[200l . AtoU 
sources have lower luminosities (~0. 001-0. 5 LEdd) than 
Z sources and have two distinct X-ray states, i.e., hard 
(energy spectra are roughly flat with power-law photon 
index near 1.7) and soft states (energy spectra follow 
exponential decrease above ~10 keV). There is also a 
"transitional" state between these two. The hard, tran- 
sitional and soft states of atoll sources are also often re- 
ferred to as the "extreme island" , "island" , and "banana" 
states/branches, respectively. Z sources only have soft 
spectra, but there are three distinct branches with dif- 
ferent spect ral and timing behav iors. Recentlv. ILin et al.l 
(|2009f l and IHoman et al.l ()2010( ) analyzed a transiently 
accreting NS, XTE J1701-462, which changed from Z- 
source behavior to atoll-source characteristics as the out- 
burst decayed from near/super Eddington luminosity 
(^Edd) to almost quiescence. These results confirmed 
that the behavior that the differences in the properties 
of the two classes are due to their different mass accretion 
rates. 

The spectral modeling of NS LMXBs has been contro- 
versial for a long time (see iB arret ,2001i . for a review) . 



The continuum of soft-state spectra in both atoll and Z 
sources are generally described by two-component mod- 
els that include a soft / therm al and a hard / Comptonized 
component (e.g., iBarret et a l. 2000: lOosterbroek et al.l 
[2OOII: iDi Salvo et al.l 120001: llaria et aLri200'5D . and there 
have been two competing models, often refer red to as 
the Eastern model (after Mitsuda ct al. 1989) and the 
Western model (after White et al. 1988), with different 
choices of the thermal and Comptonized components. In 
the Eastern model, the thermal and Comptonized com- 
ponents are described by a multicolor disk blackbody 
(MCD) and a Comptonized blackbody, respectively. In 
the Western model, on the other hand, the thermal com- 
ponent is a single-temperature blackbody (BB) from the 
boundary layer, and there is Comptonized emission from 
the disk. In the hard state, the spectra are dominated 
by a hard/ Comptonized component, but a soft /thermal 
component is generally sti l l required ([C hristian & Swan. 
199 7"; Barret ct al. '20001 iChurch fc B alucinska- Churc 
2001: Gicrlinski & Doni 120021 ). ILin et ah, (2007! ) im- 
plemented the commonly used two-component models 
for two classical transient atoll sources, i.e., Aql X-1 
and 4U 1608-52, and outlined the problem of model de- 
generacy for accreting NSs, not only from the choices 
of the thermal components, but also from the detailed 
description of Comptonized components (i.e., the scat- 
tering corona geometry, the seed photon temperature, 
etc.). The physical interpretation of the spectral evo- 
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lution of these atoll sources inferred from these differ- 
ent models also varies substantially. However, none of 
the tested models produced results similar to those of 
black- hole X-ray binaries, i.e., Lx oc tracks for the 
MCD component and weak Comptonization for the soft- 
state spectra when the integrated rms variability in the 

ower density spectrum (0.1-10 Hz) is only a few percent 

Remillard fc McClintocM[2006l) . 

ILin et al.l (I2007D devised a hybrid model for atoll 
sources, based on a detailed study two frequently re- 
curring atoll-type transients. This model uses a BB to 
describe the boundary layer plus a broken power law 
for the hard state, and two strong thermal components 
(MCD and BB) plus a constrained broken power law 
(when needed) for the soft state. This choice for the 
soft state offers a weak-Comptonization solution that dif- 
fers from the strong Comptonization solution of the two- 
component models that had previously dominated the 
literature. The results of the application of the hybrid 
model can be summarized as follows: both the MCD 
and BB evolve approximately as Lx oc T^, the spec- 
tral/timing correlations of these NSs are aligned with 
the properties of accreting black holes, and the visible 
BB emission area is very small but roughly constant 
over a wide range of Lx th at spans both the hard and 
soft states. ILin et all (120091 ) applied this X-ray spectral 
model to XTE J1701-462, and similar results were ob- 
tained for the observations when the source displayed 
atoll-source behavior. Deviations of the MCD from the 
Lx oc T"* track were observed when the source was bright 
and behaved as a Z source, with the inner disk showing 
a luminosity-dependent radius increase. This was inter- 
preted as an effect of maintaining the local Eddington 
limit in the inner disk edge as the mass accretion rate 
varies. 

The hybrid spectral model is still empirical, especially 
the modeling of Comptonization. Moreover, this model 
has only been applied thus far to extensive data ob- 
tained with the Rossi X-ray Timing Explorer (RXTE ). 
RXTE has two pointing instruments, which cover the 
energy range from ^^2.5 to 250 keV. However, there is 
emission from X-ray binaries below this energy range, as 
the characteristic temperatures of the thermal compo- 
nents are normally below 3 keV. Thus, it is important to 
test this model using broad-band spectra that extend to 
photon energies below the sensitivity range of RXTE . 

In this paper we inv estigate the persisten tly bright 
atoh source, 4U 1705-44 (jHasinger fc van der K hs 1989), 
which was observed s even times in 2006-2008 by Suzaku 
(|Mitsuda et al.l [20071 ). One of the important features 
of Suzaku is its broad energy band (0.2-600 keV). We 
also analyzed two observations of 4U 1705-44 made with 
BeppoSAX (0.1-300 keV) in 2000. Both Suzaku and 
i?eppo5'^X additionally provide better energy resolutions 
than RXTE (~ 2%, 8%, 18% at 6 keV (FWHM) respec- 
tively), and this capability can be used to better resolve 
the broad Fe emission lines. Broad Fe lines are commonly 
seen in X-ray binaries and provide another tool for inves- 
tigating the accretion flow aroun d comp act objects (e.g., 
IMilleril2007t ICackett et ani2?)08L I2009bh . 

Timing studies for 4U 1705-44, including the 
findings of kilohertz quasi-periodic oscillations (kHz 
QPOs), have b een car r ied o u t using ob servations 
with RXTE (,Ford et al.l [I998I: iBarret fc Olive .20021: 



lOlive et aLll2003[ ). Spectral studies of this source have 
also been carried out, using differ ent kinds of X-ray de- 
tectors and sp ectral models (e.g., 'Barret fc Olivel l2002l: 
|Di Salvo et"all ,2 005: Fiocchi et al. 2007; Pirain o et al l 
2007: ' HomaTet al.ll2009t IReis et al.1120091 : idi Salvo et al.l 



200_£). Several of these authors also reported the detec- 
tion of a broad relativistic Fe line from this source. In 
this paper we concentrate on the spectral properties of 
this source. We describe our data analysis in !j2l where 
we also present the long-term light curves and color-color 
diagrams. We perform detailed spectral modeling in !j3l 
for which we provide our physical interpretations in SJH 
Finally we give our summary and discussion. 

2. OBSERVATIONS AND DATA REDUCTION 

The long-term light curve of 4U 1705-44 is shown in 
Figure [TJ The gra y solid line is frorn . the RXTE All- 
Sky Monitor (ASM: iLevine et allfTOOl) . and we can see 
that the source displays persistent X-ray emission, with 
one or two strong intensity cycles per ye ar. These cy- 
cles c orrespond to state-transition cycles (jHoman et al.l 
120091) . Discrete plot symbols show the time and spectral 
states during different pointed observations. The data 
reduction and spectral analyses are described below. 

2.1. Suzaku Data 

Suzaku made seven observations of 4U 1705-44 in 2006- 
2008. Detailed information of these observations is given 
in Table [H Bo th the X-ray I maging Spectrometer (0.2- 
12 keV, XIS; Kovama et al.li2007.) and the Hard X -rav 
Detector (10-600 keV, HXD: iTakahashi et al.ir2007[) in- 
struments were used during these observations. There 
are four XIS detectors, numbered as to 3. XISO, 2 
and 3 all use front-illuminated CCDs and have very sim- 
ilar responses, while XISl uses a back-illuminated CCD. 
XIS2 was damaged in 2006 November, and its data are 
analyzed only for the first three observations. The HXD 
instrument includes both PIN diodes (10-70 keV) and 
GSO scintihators (30-600 keV). Both the PIN and GSO 
are collimated (non-imaging) instruments. 

We reprocessed each observation using the aepipeline 
tool provided by Suzaku FTOOLS version 15 and apply- 
ing the latest calibration available as of 2010 February. 
We then applied the publicly available tool aeattcor . si 
by John E. Davis to obtain a new attitude file for each ob- 
servation. This tool corrects the effects of thermal flexing 
of the Suzaku spacecraft and obtains more accurate esti- 
mate of the spacecraft attitude. For all our seven obser- 
vations, the above attitude correction produced sharper 
PSF images. With the new attitude file, we updated the 
XIS event files using the FTOOLS xiscoord program. 

Since 4U 1705-440 is a relatively bright source, win- 
dow and burst options were adopted during each ob- 
servation to limit the effects of event pile- up (Table [Ij . 
Despite this, pile- up was still present in the image cen- 
ter. We estimated the pile-up fractions at different 
positions of the CCD using the publicly available tool 
pileup_estimate . si by Michael A. Nowak. The pile- 
up fraction refers to the ratio of events lost via grade or 
energy migration to the events expected in the absence of 
pile-up. The unfiltered pile-up fractions integrated over 
the whole CCDs are about 10-15% for all observations 
except for the observation 401046010 (-3%). Table [1] 
lists the radii of the central circular regions that contain 
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Fig. 1. — The long-term RXTE /ASM light curve of 4U 1705-44. The various symbols represent different spectral states from different 
pointed observatories; BeppoSAX/ Suzaku hard state (diamond/square) and BeppoSAX/ Suzaku soft state (triangle/cross). 



TABLE 1 

The Suzaku observations of 4U 1705-44 in 2006-2008 



Observation ID 


401046010 


401046020 


401046030 


402051010 


402051020 


402051030 


402051040 


Spectral ID 


suzl 


suz2 


suz3 


suz4 


suz5 


SMZ6 


suz7 


Observation Date 


2006/08/29 


2006/09/18 


2006/10/06 


2007/09/05 


2007/10/08 


2008/02/20 


2008/03/18 


Exposure of XIS/PIN (ks)=' 


10.9/12.8 


12.7/12.1 


13.3/12.4 


2.2/8.2 


3.9/14.5 


18.4/17.3 


3.0/11.1 


XIS Detectors Analyzed 


12 3 


12 3 


12 3 


1 3 


13 


13 


1 3 


Window Option 


1/4 


1/4 


1/4 


1/4 


1/4 


1/4 


1/4 


Exposure in Burst Option (s) 


1.6 


2.0 


2.0 


0.5 


0.5 


2.0 


0.5 


XISO Count Rate (cts/s)'' 


18.1 


103.7 


45.5 


339.2 


225.9 


107.1 


229.0 


Soft Color'^ 


0.81 ±0.03 


0.72 ±0.02 


0.64 ±0.02 


0.81 ±0.03 


0.73 ±0.02 


0.66 ±0.02 


0.79 ±0.03 


Hard Color'^ 


0.34 ± 0.03 


0.11 ±0.01 


0.09 ±0.01 


0.08 ±0.01 


0.08 ±0.01 


0.09 ±0.01 


0.09 ±0.01 


Radius of Central Region 
















With >5% pile-up (pixels) 





55 


32 


52 


39 


60 


35 


Radius of Central Region 
















With >10% pile- up (pixels) 





39 


21 


35 


24 


42 


24 


Spectral state 


hard 


soft 


soft 


soft 


soft 


soft 


soft 



^ The dead time and burst clock options have been taken into account. All XIS detectors have the same exposure time, and the values 
given are for one detector only. 

^ Total time-averaged count rates, corresponding to the fully integrated PSF. 

^ The average and standard deviation of the colors for each observation based on 128s data bins. 

fractions of all annular regions are '^3% and '^5%, re- 
spectively. Using the models for the soft-state spectra, 
which include MCD and BB components, we find that 
the spectral fitting results using 10% pile-up exclusion 
regions show systematic decrease in the soft-component 
(MCD) flux and increase in the hard-component (BB) 
flux, by about 3%, compared with the results using 5% 
pile-up exclusion regions. The differences in most cases 
are within the error bars at a 90% confidence level, and 
the conclusions of this paper hold for either case. For 
simplicity and increased accuracy, we only show results 
using the 5% pile-up exclusion regions below. 

4U 1705-44 is in the direction toward the Galactic 
ridge, and the background consists of non- X-ray (par- 
ticle) background, absorbed cosmic X-ray background, 
and Galactic ridge emission. Even though our source 
can be regarded as a point source, there is no region of 
the detector plane that is free from the source emission to 
estimate the pure background. This is because the PSF 
of the XIS is quite spread out and our source during these 
seven observations was bright enough to dominate over 
the background (1-10 keV) throughout the 1/4 window. 
We compared results using two methods of estimating 
the background. In the first method, we set the back- 
ground region to be the whole 1/4 window excluding a 
circular region of radius 350" around the source. In the 
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Fig. 2. — Color-color and hardness-intensity diagrams of 
4U 1705-44 based on Suzaku observations in 2006-2008, with bin 
size 128 s. The squares are from observation 401046010, and the 
crosses from the other observations. For the definitions of the soft 
and hard colors, see the text. One-tr statistical error bars are also 
shown. 

most of the XIS CCD pixels with local pile-up fractions 
that exceed 5% and 10%, respectively. We used annu- 
lar regions to extract spectra (circular regions are used 
for observation 401046010). The outer radius was set to 
be 120 pixels, while two cases of inner radii were used, 
corresponding to the 5% and 10% pile-up exclusion re- 
gions, respectively. The corresponding integrated pile-up 
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second method, we estimated the non- X-ray background, 
which varies with time, using the xisnxbgen tool based 
on the night Earth data by Suzaku. We estimated the X- 
ray background, inchiding the cosmic X-ray background 
and Galactic ridge emission, using observation 100026030 
by Suzaku. It is specific for observing the background 
emission around the supernova remnant RX J1713-3946, 
and the pointing direction of this observation has a 184" 
offset from 4U 1705-44. Both of the methods turn out to 
give very similar spectral fit results, and we only show 
results using the second method below. 

The response files of the XIS for each observation 
were generated using the xisresp script which uses the 
xisrmf gen and xissimarf gen tools (specifying 1% ac- 
curacy). They take into account the time variation of 
the energy response and the specific extraction region 
for each observation and each XIS detector. As the re- 
sponses of XISO, 2 and 3 on the whole are very similar, 
we combined their spectra and responses using the script 
addascaspec. 

We also extracted the PIN spectra. The non X-ray and 
cosmic X-ray backgrounds are taken into account. The 
non X-ray background is calculated from the background 
event files distributed by the HXD te am. The cosm ic X- 
ray background is from the model bv iBoldtJ ()1987D . and 
its flux is about 5% of the background for PIN. The re- 
sponse files provided by the HXD team were used. The 
GSO data were not used, considering the large uncer- 
tainty in calibration and low signal to noise ratios above 
40 keV. 

The CD/HID of these observations are shown in Fig- 
ure [21 We defined soft and hard colors as the ratios of the 
count rates in the (3.6-5.0)/(2.2-3.6) keV bands and the 
(11.0-20.0)/(5.0-8.6) keV bands, respectively. The count 
rates of the lowest three energy bands were from XISO, 
1 and 3 combined. We first obtained the count rates us- 
ing the 5% pile-up exclusion regions and then converted 
to the value corresponding to the whole integrated PSF. 
They are background subtracted and deadtime in burst 
option corrected. XIS2 was not used for this because 
it was not on for all observations. The count rates in 
the energy band 11.0-20.0 keV were from the PIN with 
the background subtracted and the deadtime correction 
made. The CD and HID in Figure [U use 128-s data. 
The intensity is from the energy band 2.2-8.6 keV. Two 
type-I X-ray bursts were found in observation 401046010, 
and data around them are not included in Figure [2] or 
our spectral analysis. The data points with hard color 
larger than 0.2 (square symbols) are all from observation 
401046010, indicating that only this observation was in 
the hard state while all other observations (cross sym- 
bols) were in the soft state. 

Most observations show little variation (the average 
and standard deviation of colors are given in Table [1]) , 
and we created one spectrum each observation for spec- 
tral modeling. Each spectrum has three instrumental 
components, i.e., XIS023 (combination of XISO, 2, 3), 
XISl, and PIN. In the end, we have seven spectra from 
Suzaku observations, and they are denoted as suzl-suz7 
hereafter (Table [T]) . We rebinned the spectra by factors 
of 8 and 16 for energies below and above 2.55 keV, re- 
spectively, and further rebinning was made so that every 
bin has at least 40 counts and minimization criterion 
can be used in our spectra fitting. 



TABLE 2 

The BeppoSAX observations of 4U 1705-44 



Observation ID 


21292001 


21292002 


Spectral ID 


saxl 


sax2 


Observation Date 


2000/08/20 


2000/10/03 


Exposure of LECS/MECS/PDS (ks) 


20.6/43.5/20.2 


14.8/46.8/20.7 


MECS2 Count Rate (cts/s) 


27.9 


4.5 


Spectral state 


soft 


hard 



2.2. BeppoSAX Data 

There are two pointed observations of 4U 1705-44 
with BeppoSAX, one on 2000 August 20 in the soft 
state and the other on 2000 Oc tober 3 in the hard 
state (Table [2l IMocchi et al.|[2007[ ). The pubhcly avail- 
able data are from three narrow field instruments: the 
Low E nergy Concentrator Spectrometer (0.1-10 keV, 
LECS; iParmar et aD Il997| ). the Medium Energy Con- 
centr ator Spectrometer (1.3-10 keV, MECS: lBoella et aD 
119971 ) and the Phoswich I 3etection System (15-300 keV, 
PDS: lFrontera et al.|[l997l ). There are three MFCS units 
(MECSl, 2, 3), but no data from MECSl are available 
during these two observations. Thus we used only data 
from MECS2 and 3. We extracted two spectra, one for 
the soft-state observation and the other for the hard- 
state observation, and they are denoted as saxl and 
sax2 hereafter (Table [2). The LECS and MFCS data 
were extracted from circular regions of 8' radius centered 
on the source position. As our source is in the direc- 
tion of the Galactic ridge, we cannot use "blank fields" 
measurement for background subtraction for the LECS. 
Instead, we used th e semi-annuli method described in 
IParmar et all H19990 . For the MEGS, we used the "blank 
fields" method for the MEGS as described in the instru- 
ment analysis guide. The PDS spectra were also ex- 
tracted, with the background rejection method based on 
fixed Rise Time thresholds. The background for the PDS 
spectra was obtained using observations during off-source 
intervals. All the spectra for each instrument were finally 
rebinned using the publicly available template files to 
sample the instrument resolution with the same number 
of channels at all energies. 

3. SPECTRAL MODELING 

3.1. Spectral Models and Assumptions 

We fit all nine spectra, suzl-suz7 from Suzaku and 
saxl-sax2 from BeppoSAX. For Suzaku, we jointly fit 
spectra from XIS023, XISl and the PIN. We used the 
energy bands with good calibration and high signal to 
noise ratios for each instrument: 1.2-1.7, 1.9-2.2 and 2.3- 
10 keV for XIS detectors, and 11.0-40.0 keV for the PIN. 
Their relative normalizations were left free, and the spec- 
tral fit results quoted later are all from XIS023 (results 
from XISl differ by < 3% generally). For BeppoSAX 
spectra, we jointly fit the LECS, MEGS2, MEGS3, and 
PDS. We utilized the 1.0-3.5 keV energy band for the 
LECS, 1.7-10.0 keV for the MFCS, and 15.0-40.0 keV 
for the PDS (15.0-150.0 keV for the hard-state spectrum 
sax2). Their relative normalizations are also left free, 
but that of PDS relative to MEGS was constrained to 
the range 0.77-0.93, a 90% confidence interval advised 
by the instrument analysis guide. The spectral fit re- 
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Fig. 3. — Examples of unfolded spectra from different states using 
different models. The total model fit is shown as a black solid line. 
The MCD component (if included) is shown by a red dotted line, 
the BB component by a blue dashed line, the PL/CPL/SIMPL 
component by a green dot-dashed line, and the Fe line (modeled 
by the diskline model) by a cyan triple-dot-dashed line. For the 
SIMPL(MCD)-|-BB model, the MCD component shown is the un- 
scattered part, and the plotted SIMPL component refers to the 
scattered part. The hard-state spectrum is from sax2, and the 
soft-state spectrum from suz2. 

suits quoted later are all referenced to MECS2 (results 
from MECS3 differ by < 2%). The fit of the Crab Neb- 
ula from the MECS using a single power law gives the 
photon index 2.1 and normalization 9.23, while the XIS 
gives 2.1 and 9.55, respectively. Thus the normalizations 
from both observatories appear to differ by less than 5%. 
For all spectral fits, we set the model systematic error to 
be 1%. 

For the soft-state spectra, we tested three models for 
the continuum spectra: MCD+BB, MCD-^BB-hPL, and, 
SIMPL(MCD)-fBB, respectively, where PL is a power 
law and SIMPL is a simple Comptonization model by 
iSteiner et all pOOl ). MCD (diskbb in XSPEC) has two 
parameters: the temperature fcTMCD at the apparent in- 
ner disk radius i?MCD , and the other is the normalization 
^MCD = (-RMCD,km/£'iokpc)^cosi, where -Diokpc is the 
distance to the source in unit of 10 kpc and i is the 
disk inclination. BB (bbodyrad in XSPEC) assumes an 
isotropic blackbody spherical surface with radius Rbb 
and has two parameters, i.e., the temperature kT^s and 



the normalization Nqq = (i?BB,km/-Diokpc)^- SIMPL (in 
XSPEC12) is an empirical convolution model of Comp- 
tonization in which a fraction of the photons from an 
input seed spectrum is scattered into a power-law compo- 
nent. This model has only two free parameters, i.e., the 
photon power-law index Fsimpl and the scattered frac- 
tion fsc- In addition, there is a flag parameter to con- 
trol whether all the scattered photons are up-scattered 
in energy or are both up- and down-scattered. We spec- 
ified that all the scattered photons are up-scattered in 
energy. We found that inclusion of down-scattering only 
changed our results within error bars. We assumed that 
the Comptonization seed photons are from the disk. The 
fit results of the MCD and BB are consistent within error 
bars if we assumed both the MCD and BB contribute to 
the Comptonization seed photons instead, thanks to low 
Comptonization for all our soft-state observations. The 
best-fitting photon power-law index Fsimpl turns out to 
be high (> 4) in most cases, a reg ime where the model is 
not suitable (jSteiner et al.l 120091) . Thus we constrained 
LsiMPL to be less than 2.5, a value typically seen in the 
black-hole cases. No constraint on the photon index FpL 
in the PL model was used. 

The two hard-state continuum spectra sax2 and suzl 
were fit with a BB plus a Comptonized component. 
We tested three choices of the Comptonization compo- 
nent: a broken power law (BPL; bknpower in XSPEC), 
a cut-off power law (CPL; cutoffpl in XSPEC), and the 
Comptonization model by rTitarchukI (j 19941) (CompTT 
in XSPEC). CompTT is an analytic model describing 
Comptonization of soft photons in a hot plasma. 

All models included an absorption component (we use 
model wabs in XSPEC). There are strong b road Fe 
lines in most spectra (see also Reis et al.' '2009'). They 
were modeled by the diskline model (Fabian et ah . 1989 ), 
which describes line emission from a relativistic accre- 
tion disk. Its parameters are: the line energy E\[^q in 
unit of keV, the power law dependence of emissivity (/3), 
the disk inner and outer radii in units of GM/c^, the disk 
inclination z, and the normalization (photons cm~^ s~^). 
Figure [3] shows some examples of unfolded spectra at 
different states using different models, with the Fe lines 
modeled by the diskline model. The top panel shows the 
hard-state spectrum sax2 using model CPL-t-BB. The 
lower two panels show the soft-state spectrum suz2 using 
models MCD+BB+PL and SIMPL(MCD)+BB, respec- 
tively. 

We scaled the luminosity a nd radius related quantitie s 
usmg a distance of 7.4 kpc (jHaberl fc TitarchukI [19951) . 
unless indicated otherwise. The flux and its error bars 
are all calculated for an energy band of 0.001-200 keV 
(1.5-200 keV for CPL/PL components) using the cflux 
model in XSPEC12. As there is little emission of the 
thermal components outside the above energy band, the 
values are essentially bolometric for thermal components. 
For the MCD component, we assume the inclination to 
be 24°, from the fitting of Fe lines (see below). 

3.2. Spectral Fit of the Soft-state Spectra and Results 

We test three models for the soft-state con- 
tinuum spectra: MCD-hBB, MCD+BB+PL, and 
SIMPL(MCD)+BB. Our practice is to use a fixed value 
of the interstellar column density (Nh) for all spectral 



6 



Lin et al. 



TABLE 3 

Spectral modeling results of soft-state observations using MCD+BB+PL+diskline 



fcTMCD (keV) Wmcd feTBB (keV) ATbb Tpl Wpl -Bii„o (kcV) EW (eV) xlH ix.Edd 

suz3 0.82 ±0.02 114 ± 13 1.87 ± 0.03 6.43 ± 0.53 2.71 ± 0.080.17 ± 0.02 6.87_o.i2 -4.55 ± 0.50 160 ± 24 1.27(703) 0.044 ± 0.001 

suz6 1.07 ±0.02 101 ± 6 2.12 ± 0.04 8.05 ± 0.57 2.77 ± 0.070.33 ± 0.02 6. 92_o og -4.06+°;!^ 176^23 1.05(708) 0.098 ± 0.002 

suz2 1.17 ± 0.03 63 ± 4 2.09 ± 0.04 9.53 ± 0.91 2.65 ± 0.070.30 ± 0.02 6. 94_o og -3.64 ± 0.27 160 ± 20 1.04(708) 0.099 ± 0.002 

saxl 1.31 ±0.03 68 ±5 2.08151 14.46t° jg 2.69t°;o5 0.62 ± 0.03 6. 97_o o5 -3.59 ± 0.20 156^" 1.08(217) 0.163 ± 0.003 

suz5 1.37 ±0.02 88 ± 5 2.25 ± 0.03 11.73 ± 0.892.72 ± 0.090.38 ± 0.04 6. 95_o oei -3.45 ± 0.28 154tll 1.00(706) 0.192 ± 0.004 

suz7 1.57 ±0.06 48 ±5 2.3 ± 0.04 11.5 ± 1.47 2.67 ± 0.080.44 ± 0.04 6. 89_o.os -S.SStoH 165 ± 260.99(706) 0.198 ± 0.007 

suz4 1.65 ± 0.07 61^6 2.33 ± 0.05 15.51 ± 2.48 2. 77t° Jj 0.50 ± 0.06 6. 97_o.o5 -3.4 ± 0.30 list^^ 1.00(702) 0.285 ± 0.012 



Note. — They arc listed in order of the souree total luminosity. Sec ti3. ll for the meaning of each parameter. The normalizations 
of the MCD and BB models are based on the assumption of the distance to the source of 10 kpe. For the diskline model, the 
inclination is fixed at 24° , and inner disk radius at 6GM/c^ (see text). EW is the equivalent width of the Fe line modeled by the 
diskline model, but it is not a parameter of the model. The last column is the total luminosity in unit of the Eddington luminosity 
( j3T2t . and the error bars are calculated based on simple error propagation from individual spectral components. 



TABLE 4 

Spectral modeling results of soft-state observations using SIMPL(MCD)+BB+diskline 



feTMCD (keV) AfMCD fcTBB (keV) Nbb Lsimfl fsc Biino (keV) /3 EW (eV) xlii^) ix.Edd 

suz3 0.81 ±0.01 139 ± 6 1. 86 ± 0.02 6. 67 ± 0.37 2. 5_o 15 0.06 ± 0.02 6. 86_o 12 -'i-57to,t9 163 ± 23 1.28(703) 0.044 ± 0.001 

suz6 1.04±0.01 127 ± 42.10 ± 0.03 8.63 ± 0.59 2.5_o ig 0.04 ± 0.01 6.95_o 09 -4.28 ± 0.43 173 ± 21 1.06(708) 0.097 ± 0.002 

suz2 1.10± 0.02 90±4 2.04±0.0410.82±0.812.5_o 15 0.06^° 02 6.95_oo6 -3-79 ± 0.25 165 ± 20 1.05(708) 0.100 ± 0.002 

saxl Lnt^M 12ltg 2.03 ±0.02 n-ieto.ll 2.5_o 1 0.05 ± 0.01 6.97_o 03 -3.87 ± 0.19 187^15 1-31(217) 0.165 ± 0.002 

suz5 1.39 ±0.02 88 ± 3 2.27 ± 0.03 11.19 ± 0.77 2. 5_o 89 0.02 ± 0.01 6. 93_o 06 -3-30to:il 161 ± 29 1.02(706) 0.185 ± 0.004 

suz7 1.49 ±0.03 62 ± 3 2.27 ± 0.0412.90 ± 1.18 2. 5_o 45 0.03 ± 0.01 6.9_o Qg -^■'^^-oil 169^2^ 0-99(706) 0.197 ± 0.005 

suz4 1.66 ± 0.05 65 ± 5 2.34 ± 0.05 14.94 ± 2.062.11_o 61O.OI ± 0.01 6.97_o 05 -3.37 ± 0.28 117^15 1.00(702) 0.277 ± 0.010 



Note. — Same as Table[3] but for model SIMPL(MCD)+BB+disklinc. 



fits of each model. To determine for the appropriate 
value, we fit all soft-state spectra {suz2 — suz7 and saxl) 
simultaneously with each model, while tying their iVjj 
to a common value. The best-fitting values of iVn for 
MCD+BB, MCD+BB+PL, and SIMPL(MCD)-|-BB are 
(1.69±0.01), (1.88±0.06), and (1.71±0.01) x 10^2 cm-^, 
respectively. The above values are obtained from fits 
with the Fe line region 4.0-8.0 keV excluded, although 
the fits with the Fe line modeled by the diskline model 
give very similar values. We see that models MCD+BB 
and SIMPL(MCD)+BB yield similar values of Nu, but 
they are smaller than that from m odel MCD-I-BB-I-PL . 
This is consistent with the results of lSteiner et al.l (|2009l ). 
i.e., fits with PL tend to infer higher values of A'h than 
fits with SIMPL, a systematic effect of different ways 
of handling the Comptonization. Considering this issue, 
we do not require a common value of A'h for all models, 
but use A^H = 1-71 x lO^^ cm'^ for models MCD-I-BB 
and SIMPL(MCD)-|-BB, and 1.88 x 10^2 cm'^ for model 
MCD+BB+PL. In the latter case, the choice of either 
value of Nh does not affect the spectral parameters for 
the MCD and BB components very much, but there are 
obvious differences for the PL component, as will be dis- 
cussed below. 

The spectral fit results of the thermal components 
(MCD and BB) in the soft state from all tested 
models with the Fe line included are shown in Fig- 
ure 131 and results of all spectral components are 
tabulated in Tables |3] (MCD+BB+PL+diskhne) and 
H (SIMPL(MCD)+BB+diskline). The details of fit- 
ting the Fe line with the diskline model are given 
in §3.41 The red crosses and the black triangles 



in Figure |4] are for soft-state spectra suz2-suz7 and 
saxl, respectively, and the panels from the top to the 
bottom correspond to continuum models MCD+BB, 
MCD+BB+PL, and SIMPL(MCD)+BB, respectively. 
For model SIMPL(MCD)+BB, the MCD component 
shown is the original value before scattering. The dot- 
ted lines in Figure |4] correspond to the NS burst radius 
of i?burst ~ 7.4 km (at a distance of 7.4 kpc), assuming 
Lx — ■inR^aT'^. The NS radius was derived from spec- 
tral fitting to Type I X-ray bursts of this source using 
RXTE data (see also iGottwald et aLllTOSl ). The dashed 
lines correspond to i? = 2.2 km, which is about the av- 
erage visible BB emission size. 

From Figure 21 we see that the MCD and BB compo- 
nents in the soft state roughly follow the L oc T"* tracks 
for all models, which implies relatively constant appar- 
ent emission areas. We discuss the extent of deviation 
in [jH The inner disk radius is comparable with the NS 
radius, while the visible BB emission area is about 1/11 
of the NS surface. The kTucD has values from ~0.8 to 
1.7 keV, and fcTeB from -1.8 to 2.4 keV. These values 
are roughly similar to those seen in Aql X-1, 4U 1608- 
52, and XT E J1701-462 in the atoll soft state ()Lin et al.l 
120071 12009f ). Fits with the Fe line region excluded from 
the fit generally give consistent results for the MCD and 
BB components, to within 10% for the normalization pa- 
rameter A'mcd and within 20% for the Nbb- Spectrum 
saxl gives larger differences (— 40% for A'bb), probably 
because of its narrower energy band and lower energy 
resolution (only five channels above 10 keV and none in 
the energy band of 10-15 keV). 

Results of the thermal components from both mod- 
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Fig. 4. — The luminosity of the thermal components versus their 
characteristic temperatures. Each model includes the diskline com- 
ponent to fit the Fe line. The symbols in each panel denote dif- 
ferent spectral states and different observatories: BeppoSAX hard 
(black diamond) and soft (black triangle) states, and Suzaku hard 
(blue square) and soft (red crosses) states. For the case of Model 
SIMPL(MCD)-|-BB, the MCD component shows the original seed 
spectrum luminosity (i.e., before scattering). The dotted lines cor- 
respond to the NS burst radius of 7.4 km (see the text), and the 
dashed lines correspond to Rbb = 2.2 km. Error bars at a 90% 
confidence level are also shown, but mostly are smaller than the 
symbol size. 

els MCD-I-BB-PPL and SIMPL(MCD)+BB are generally 
similar, but differences of > 20% can occur in some cases 
(e.g., 30% for A^mcd from the spectrum suz2; compare 
Tables [3H11) • Due to weak Comptonization in the soft 
state, model MCD-f BB in general also gives similar re- 
sults of the thermal components (e.g., flux differs < 5% 
and A^BB by < 20%) and acceptable reduced values 
(< 2.0; but 2.6 for saxl). The largest differences in 
best-fitting spectral parameters from model MCD-I-BB 
compared with those from models MCD-I-BB-I-PL and 
SIMPL(MCD)-fBB are from BeppoSAX spectrum saxl 
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Fig. 5. — The fit residuals in terms of sigmas with error bar values 
set to be one for different Comptonization models in combination 
with BB model, using the hard-state spectrum sax2. The red, 
black, green points are from LECS, MECS2, and PDS, respectively. 
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Fig. 6. — The energy fraction of Comptonized luminosity ver- 
sus the total luminosity, from model SIMPL(MCD)-|-BB (the soft 
state) and model CPL-f BB (the hard state). See Figure[4]for mean- 
ings of symbols. Error bars at a 90% confidence level are shown 
for the fraction of Comptonization Lx, and those for Lx/LEdd are 
not shown, but all are very small (Tables [3}0. 

{Nbb differs by -50%). 

3.3. Spectral Fit of the Hard-state Spectra and Results 

We first investigate the performance of different Comp- 
tonization models. Using spectrum sax2, which ex- 
tends to 150 keV (instead of 40 keV for spectrum 
suzl), we obtained values of x^/^^ 190.2/226, 287.1/224, 
and 339.9/225 for models CPL-PBB, BPL+BB, and 
CompTT-|-BB, respectively {Nn was allowed to have 
different values for different models). The result for 
CompTT-t-BB quoted above is just one possible solution 
with a local minimum. The best-fitting input seed 
photo n temperature rp is <0.2 keV (the cold seed photon 
model: iLin et al.ll2007[ ). All the above models give similar 
results in terms of fcTeB (~ 1-5 keV), TVbb (^ 20), and 
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TABLE 5 

Spectral modeling results of hard-state observations using BB+CPL+diskline 



Nn (10^^ cm-^)TBB (koV) 



Nb 



TcPL 



Afcp 



En 



(kcV) 



EW (oV) xli'^) 



suzl 
sax2 



suzl 
sax2 



1.71f 
1.71f 



1.88f 
1.88f 



0.94 ± 0.0217.62 ± 1.581.19 ± 0.020.056 ± 0.0026.60 ± 0.04-2.45 ± 0.28 92 ± 27 1.10(676)0.055 ± 0.001 



1.30 ± 0.03 6.06 ± 0.61 1.37 ± 0.010.119 ± 0.003 6.80_ 



-2.46 ± 0.55 93 ± 36 1.19(228)0.070 ± 0.001 



0.94 ± 0.0313.12 ± 1.621.33 ± 0.010.076 ± 0.0026.59 ± 0.04-2.30 ± 0.33 74 ± 28 1.06(676)0.049 ± 0.001 



1.24 ± 0.04 5.52 ± 0.77 1.40 ± 0.010.135 ± 0.003 6.87_ 



-2.86 ± 0.55 126 ± 370.93(228)0.070 ± 0.001 



Note. — Same as Table |3] but for model BB-t-CPL-t-diskline and for the hard-state data. The cutoff energy of the CPL is 
fixed at 44 keV. 



/ 



10"^ 0.1 1 

BB (10^" erg/s) 

Fig. 7. — The luminosity of the MCD component plus Comp- 
tonization versus the BB luminosity. The results for the soft state 
are from model SIMPL(MCD)-t-BB-f diskline, and the Comptoniza- 
tion refers to the scattered disk emission, modeled by SIMPL. 
Model MCD-f-BB-f PL-|-diskline gives similar results. The dashed 
and dotted lines mark the ratios of 1 and 0.05 of the BB luminos- 
ity versus the MCD component plus Comptonization luminosity, 
respectively. See Figure |2] for meanings of symbols. Error bars at 
a 90% confidence level are also shovirn, but mostly are smaller than 
the symbol size. 

the Comptonization fractions (> 90%). There is another 
sohition for CompTT-l-BB, which has xV'^ = 252.5/225. 
It has tq = 0.9 ± 0.1 keV (the hot seed photon model; 
iLin et aP l2007| l. For this case, the inferred BB com- 
ponent is quite different from the models above, with 
fcTBB = 0.5 ± 0.1 keV and A^bb = 256 ± 110. The fit 
residuals of all the above models are shown in Figure [5l 
We see that CPL-f BB gives the best fit over the whole 
energy range, while other models tend to have large resid- 
uals at high energy (> 15 keV). Hereafter we focus on 
model CPL+BB only. 

The fit of the CPL model to the BeppoSAX hard-state 
15.0-150.0 keV spectrum (the BB contributes little at 
energies above 15.0 keV) gives a cutoff energy of 44 ± 4 
keV. The cutoff energy cannot be well constrained in the 
case of spectrum suzl, for which we use energies only up 
to 40 keV. Thus we fixed the cutoff energy to be 44 keV 
when we fit spectra sax2 and suzl. The simultaneous 
fit of sax2 and suzl gives a best-fitting value of A'h = 
(1.89 ± 0.04) X 10^2 cm-2 from model CPL-t-BB, close 
to the value of 1.88 x 10^^ cm^^ obtained from model 
MCD-I-BB-I-PL for the soft state. For comparison, we fit 
the hard-state spectra using both Nn = 1.88 and 1.71 x 



lO'^^ cm^^, and the results are given in Table [5) We 
see no large different between these two sets of solutions. 
Using A^H — 1.88 X 10^^ cm~^ tends to give a smaller 
BB emission area and a higher photon index of the CPL 
model. The values of the photon index are more similar 
between spectra sax2 and suzl than using A'h — 1.71 x 
10^^ cm^^. Both values of A^h give high fractions of 
Comptonization (>90%). 

The results using A'h — 1-88 x 10^^ cm^^ are shown 
in Figure |4l the blue squares and the black diamonds for 
suzl and sax2, respectively. They are repeated in the 
middle and bottom panels. They are not shown in the 
top panel, because this panel is reserved for solutions free 
of Comptonization and such a model is unacceptable for 
observations of the hard state. From Figure HI we see 
that the boundary layer emission areas in the hard state 
are comparable to those in the soft state. Although L-qb 
changes over 50 times, the values of A'bb are always in 
the range of ^^5-15. 

One of our main goals of using the broad-band spectra 
is to search for detectable thermal emission from the ac- 
cretion disk in the hard state. We added the MCD com- 
ponent while analyzing the two hard-state spectra {suzl 
and sax2), i.e., using model CPL-I-MCD-I-BB, with ei- 
ther the Fe modeled with the diskline model or with the 
Fe region excluded in the fit. A^h is either fixed at 1.71 or 
1.88 xlO^^ cm-^ or left free in the fit. The best-fitting 
disk temperature tends to go below 0.2 keV, with an 
upper limit < 0.25 at a 90%-confidence level. The nor- 
malizations of the MCD component A'mcd are not well 
constrained. For all cases, the unscattered flux of the 
MCD component is < 2% of the total flux (absorbed or 
unabsorbed; 1-200 keV). Thus, if we assume a physically 
visible disk, then we cannot exclude either possibility, 
i.e., that the disk in the hard state might be truncated 
at a very large radius and/or the temperature is below 
0.2 keV. Alternatively, the disk may be rendered invisible 
by very high Comptonization. 

Figure IHl shows the fraction of Comptonized luminos- 
ity versus the total luminosity, using SIMPL (MCD) -fBB 
for the soft state and CPL+BB for the hard state. The 
total luminosity is normalized by the Eddington lumi- 
nosity LEdd, which is derived from the type I X-ray 
bursts showing photospheric radial expansion and cor- 
responds to a n average peak flux of about 4 x 10~^ 
erg cm" -2 (|Gallowav et alj I2OO60 . Figure [6] shows 
that Comptonization only constitutes < 15% of the emis- 
sion in the soft state, but > 90% in the hard state. 
For the soft-state data, the fractional contribution of 
Comptonized luminosity decreases with luminosity on 
the whole, which is consistent with the behavior of atoll- 
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type transients (|Lin et alj [20071 ) . The MCD+BB+PL 
model also gives low-Comptonization solutions for the 
soft state (< 15%), with either 1.71 or 1.88 xlO^^ 
cm"^. One main difference between these two choices of 
A^H is that fits with = 1.88 x 10^^ cm'^ give higher 
values of Tpl 2.7) than fits with Nu = 1.71 x 10^^ 
cm-2 (FpL 2.2). 

The BB component, which describes the emission from 
the boundary layer in our model, is present in both the 
hard and soft states. We can compare this component 
with the other spectral components to investigate the im- 
pact of different accretion processes in different states. 
We plot in Figure [7] the luminosity of the MCD com- 
ponent plus Comptonization (SIMPL/CPL) versus the 
BB luminosity from model SIMPL(MCD)+BB. Model 
MCD+BB+PL gives similar results. The hard-state data 
(diamond and square symbols) are from model CPL+BB. 
The dashed and dotted lines correspond to the ratios of 
1 and 0.05 of the BB luminosity versus the MCD compo- 
nent plus Comptonization luminosity, respectively. They 
are about the average values for the soft- and hard-state 
observations, respectively. If we assume that the Comp- 
tonization emission is not from the boundary layer, the 
above result means that there is a much lower portion of 
energy seen in the visible portion of the boundary layer 
in the hard state than in the soft state. Similar results 
were suggested to be a possible consequ ence of a strong 
jet in the hard state in lLin et al.l (|2007[ ). 

3.4. Relativistic Fe lines 

We first fit the Fe line with a Gaussian line. A broad 
Gaussian line with the line width a around 0.6 is gen- 
erally required, while fits with a narrow line (cr < 0.1 
keV) are mostly unacceptable, with values increased 
by > 100 for ^700 degrees of freedom in the soft state. 
The diskline model can in general improve the fits fur- 
ther, compared with a broad Gaussian line, with the 
values decreased by about 30 on average in the soft state. 
In the hard state, the diskline model still gives the best 
fits, though the improvement is much less. We also fit 
the spectra with a smeared edge (smedge in XSPEC) in 
our continuum model, with the lower limit of the thresh- 
old energy set to be the neutral Fe K edge 7.1 keV. We 
obtain fits with values larger than those using the 
diskline model by about 30 on average, and the thresh- 
old energy in most fits reaches the lower limit 7.1 keV. 
Thus the broad line feature in the spectra is probably 
not completely due to the Fe absorption edge. The more 
complicated model combining both a Gaussian line and 
a smeared edge can give fits with similar values as 
those using the diskline model, but we choose to focus 
on the results using the simpler model, i.e., the diskline 
model. 

While fitting the Fe line with the diskline models, we 
initially fit all soft-state spectra simultaneously, with the 
inclination parameter i tied to a common value. Dif- 
ferent continuum models turn out to give quite simi- 
lar best-fitting values of i: 24.3 ± 0.8°, 23.3 ± 1.0°, and 
24.2 + 0.8° from models MCD+BB, MCD+BB+PL, and 
SIMPL(MCD)+BB, respectively. Thus we fix the disk 
inclination to be i = 24°. For most of the soft-state 
spectra, the inner disk radius inferred from the diskline 
model reaches 6 GAI/c^, the innermost stable circular 
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Fig. 8. — Fe lines of all spectra, fit by the diskline model. The 
upper two panels are for the hard state (spectra suzl and sax2), 
while others for the soft state. The soft-state continuum is fit 
by model SIMPL(MCD)+BB, while model MCD+BB+PL gives 
very similar results. From top to bottom, the source luminosity 
increases, except for the hard-state spectra suzl and sax2, which 
are put on the top panels and have luminosities between that of 
suz3 and suz6. 

orbit (ISCO). Thus in the final fitting, we fix it to be 
6 GM j (? . For the hard-state spectra, the Fe lines are 
much weaker, and the best-fitting inner disk radius in 
the diskline model is quite uncertain. The lower error 
bar reaches 6 GMjt? for the two observatories using both 
A^jj = 1.88 and 1.71 x 10^^ cm~^. For simphcity, we then 
fixed the inner disk radius in the diskline model to be 
6 GM/c^ for the hard-state spectra to derive the final 
results. The obtained Fe line fiux and equivalent width 
are not sensitive to such details. 
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Fig. 9. — Dependence of Fe line flux on the BB flux (upper panel) 
and Comptonizatlon flux (SIMPL/CPL components; lower panel). 
The solid lines mark the ratio of 4% of the Fe line flux over the BB 
flux or the Comptonizatlon flux. The Fe line flux is bolometrie, but 
the BB flux and Comptonizatlon flux are integrated over 6.4-200 
keV. See Figure |4] for meanings of symbols. Error bars at a 90% 
confidence level are also shown, but mostly are smaller than the 
symbol size. 

The results of fitting Fe lines with the diskline model 
are shown in Figure HI It can be seen that the Fe emission 
line is detected in all spe ctra, and most of them show 
a broad feature (see also iReis et aD 120091 ) . Tables [SHS] 
contain best-fitting parameters of the Fe line for each 
spectrum, with the equivalent width also included. The 
best-fitting value of iSune reaches the allowed upper limit 
6.97 keV in most cases. The equivalent width is around 
170 eV for most of the soft-state spectra, and it is around 
100 eV in the hard state. The power-law dependence of 
emissivity (3 is generally steeper in the soft state (~ —3.5) 
than in the hard state —2.5). 

Fe emission lines in X-ray binaries are the most obvious 
signature of an accretion disk irradiated by an external 
source of hard X-rays, due to a combination of high fluo- 
rescent yield and large cosmic abundance (pillcr 2007j). 
To investigate the irradiation source of the Fe lines, we 
show in Figure IHl the dependence of the bolometrie Fe 
energy line flux on the BB energy flux (upper panel) 
and Comptonizatlon (SIMPL/CPL) flux (lower panel), 
both integrated over 6.4-200 keV. The soft-state data 
show that the Fe line flux increases monotonically with 
the BB flux, but has no clear dependence on the Comp- 
tonizatlon flux. This might imply that it is the bound- 
ary layer emission that illuminates the accretion disk and 
p roduces the Fe Hne, in agreement with the conclusions 
of lCackett et al.l (|2009aD . The Fe line energy flux is about 
4% of the BB energy flux (solid lines). 



The hard-state data in Figure[5] (diamond/square sym- 
bols) in the upper panel show that there is much higher 
Fe line flux relative to the BB continuum in the hard state 
compared to the soft state. This might indicate that the 
BB emission is not the (major) source illuminating the 
accretion disk to produce the Fe line in the hard state. 
From the lower panel, we can see that there is strong 
Comptonizatlon flux available in the hard state. How- 
ever, using the reference lines and comparing the upper 
and lower panels, we see that the Comptonizatlon emis- 
sion is not as efficient as the BB emission to illuminate 
the accretion disk and produce the Fe line. This can 
be explained if the Comptonizatlon emission is located 
farther away from the accretion disk than the bound- 
ary layer. We also caution that some BB flux may be 
screened from the observer's view by the accretion disk. 

4. DISCUSSION 

4.1. Global View of Two Spectral States 

We have shown that 4U 1705-44 exhibits two distinct 
spectral states. The hard-state spectra span from '^1 
to 150 keV, while the soft-st ate spectra are most ly con- 
flned below 40 keV (see also iFiocchi et al.ll2007l ). Such 
two distinct s tates have be e n observed in other atoll- type 
NS LMXBs (jBarretl [200l Ivan der Klid [20061) . On the 
whole, these two spectral states are very similar to the 
thermal and hard states (but not the steep power law) 
of accreting black holes ([Remillard fc McClintockl [20061 : 
iDone et all 120071 ). In addition, atoll- type NS LMXBs 
and accreting black holes also have many similarities in 
timing properties. All this motivates the speculation 
that the accretion process is quite similar b etween atoll- 
type NS LMXBs and accreting b lack holes ()Barretll200ll: 
IDone et al.ll2007t ILin et al.ll2007D . This implies that the 
presence of the boundary layer in atoll-type NS LMXBs 
probably does not strongly interfere with the accretion 
disk, at least at luminosities below ^0.5 L-^dd- How- 
ever, the boundary layer is still an important spectral 
component in atoll-type NS LMXBs. Figure [7] shows 
that the boundary layer emission flux is comparable with 
that from the disk. Because the bou ndary layer is hotter 
than the disk ([Mitsuda et al.|[T984D . the soft-state spec- 
tra tend to cover a broader energy band and extend to 
higher energies in atoll-type NS LMXBs than in the ther- 
mal state of accr eting black holes (compar e our F igure B] 
and Figure 2 in iRemillard fc McClintoci] (|2006( )). The 
relatively weak contribution of the boundary layer in the 
hard state, as shown in Figure [71 can be inte rpreted as 
the c onsequence of mass ejection in this state ()Lin et al.l 
[200l . This seems to be in alignment with the jet model 
for the hard stat e in accreting black holes ([Fended 120061 : 
IDone et al.l[2007[ ). 

4.2. Inner Disk Radius 

In ILin et al.l (|2007f ). our spectral modeling of two atoll 
sources Aql X-1 and 4U 1608-52, using the MCD+BB 
model plus weak Comptonizatlon showed the MCD be- 
havior close to the Xmcd oc Tmcd track. In contrast, 
classical two-component models resulted in constant or 
slightly decreasing temperatures of the thermal compo- 
nent (MCD or BB) with increasing luminosity. We note, 
however, that in that study we used RXTE data and the 
MCD component in the low-luminosity soft state could 
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not be well constrained. In this study, we used Suzaku 
and BeppoSAX spectra, which extend the energy range 
down to 1 keV, and the MCD components are well con- 
strained for all soft-state spectra. As the result of the 
better low-energy coverage, we did not need to put addi- 
ti onal constraints o n the PL component, as was necessary 
in iLin et al.l (|2Q07f ). This motivates us to evaluate more 
quantitatively how closely the MCD component follows 
the Lmcd oc T^cd track. 

Our spectral modeling of broad-band spectra of 4U 
1705-44 shows a track that is flatter than the Lmcd oc 
-^MCD relation, with nearly an order of magnitude varia- 
tion in luminosity. There is a slight decrease of -Rmcd at 
higher luminosity. Spectral fitting with Comptonization 
modeled by PL suggests Lmcd oc Tmcd"^' while fitting 
by SIMPL results in Lmcd oc Tmcd"^- Such deviations 
have been seen in several black-hole X-ray binaries, and 
it is normally believed to be due to a spectral hardening 
effect, instead of a real change in the inner disk radius 
jShafee et al."200^ iDavis et all 120061: iMcClintock et al.l 
[2007, 2009). Spectral hardening arises when the electron 
scattering dominates over absorption as the main opacity 
source. In such a situation, the local specific flux in the 
disk appears as a simple dilute blackbody with a color 
temperature high er than the effective temperature by a 
factor of /col (Shi mura fc Ta kahara 1995). This factor 
slightly increases with luminosity/temperature. Based 
on a simple analytic estimate of the hardening factor, 
IDavis etall H2006I ) suggested a Lmcd oc Tj^cd relation. 
It should be noted that the extent to which the hardening 
factor depends on luminosity can vary with the inclina- 
tion, the mass of the compact object, etc. A detailed 
numeric simulation to obtain how the hardening factor 
behaves for an accreting NS or direct accretion model- 
ing to infer the real inner disk radius with more realistic 
spectral models incorporating many effects could provide 
further insights. We further note that there are possibly 
other factors causing the above deviation. We cannot 
exclude the possibility of real change of the inner disk 
radius in some part of the soft state. It is also possible 
that the above deviation is due to our simple descriptions 
of the weak Comptonization and the boundary layer. 

4.3. Constraint on the Magnetic Field in 4U 1705-44 

We perform a rough estimate of the magnetic field in 
4U 1705-44, under the assumption that the disk is trun- 
cated at the ISCO in our soft-state observations. This 
requires the magnetic field to be dynamically unimpor- 
tant for these soft-state observations. That is, the Alfven 
radius ta, the radius at which the magnetic pressure is 
roughly the sum of the ram and gas pressure, should be 
smaller than the IS CO. Based on Equations 6.19-6.20 in 
iFrank et afl ()1985D . we have 

{jW^y'^'iT^f' km, (1) 

where B is the magnetic field strength at the surface 
of the NS and /ca is the correction from the spherical 
accretion to disk accretion and is about 0.5. This for- 
mula assumes a dipole magnetic field. We further assume 
A/ns = 1.4 Mq (such that i?in = 12.4 km), and Rm = 10 



km. The above expression shows that the obtained value 
of B will only weakly depend on the NS mass (a power 
of 1/4). Using L = lO^^ erg/s, fr om the faintest soft- 
state spectrum suz3, and the constraint ta < -Rin, we 
find S < 1.9 X 10^ G. We can also assume that ta is less 
than i?NSj which would decrease the above limit by 30%. 

4.4. The Boundary Layer 

We see that the apparent emission of the boundary 
layer, modeled by BB, roughly follows Lbb oc Lgg, from 
the hard to the soft states with the luminosity of the 
boundary layer covering ~0. 003-0. 12 LEdd- The appar- 
ent area of the boundary layer is about 1/11 of the the NS 
surface, as inferred from Type I X-ray bursts. The actual 
area of the boundary layer might be larger due to a spe- 
cial geometry of the boundary la yer, which has b een at- 
tributed to be an equatorial be lt (|Lin et al.ll2007[ ). Using 
Equation 3 in ILin et aTl (j2007l ) and using the inclination 
obtained from the Fe line fit (z = 24°), the latitude range 
(from the NS equator) of the boundary layer is about 
13°, corresponding to an emission area of 23% of the 
NS surface. Accordingly Lbb should be 1.6 times larger 
than shown in § [3l which assumes isotropic emission of 
the boundary layer. Thus, LBB/LMCD-i-Comptonization is 
not around 1 in the soft state, as shown in Figure [71 
but is about 2.6. If ? = 60°, the latitude range of the 
boundary layer is about 7°, corresponding to an emis- 
sion area of 12% of the NS surface. In this case Lbb 
should be 30% larger than shown in § [3l Lmcd should 
increase by 80% (i.e., adjusting for i — 60°) so that 
LBB/LMCD+Comptonization would be arouud 0.7 in the soft 
state. 

Whether the above results imply that the real bound- 
ary layer area is small and nearly constant further de- 
pends on the radiative transfer process in the atmosphere 
above the boundary layer, i.e., the hardening effect as 
discussed above for the disk spectra. The small BB emis- 
sion area is reminiscent of the well known spectral mod- 
cling problem of the NS thermal emission in quiescence, 
i.e., the BB fit of its thermal component produced in- 
ferred radii too small for theoretical NS size estimate, 
whereas models taking into account the radiative transfer 
in the hydrogen atmosphere give radius estimate much 
closer to theoretical expectation of the size of NSs (e.g., 
iRutledge et"anil999[ ). However, all our observations are 
quite bright (> 10'^'' erg/s), and the emission should be 
due to active accretion. At such a high acc retion rate, a 
pure hydrogen atmosphere is not expected ([Brown et al.l 
Il998f ). and the above problem might not apply to our 
case. Our conclusion of the small size of the boundary 
layer is based on the assumption that the modification 
of the bursting atmosphere on burst emergent spectra 
is similar to that of the boundary layer emission. This 
assumption might be valid if most of the heat in the 
boundary layer is generated in a layer as deep as that 
for burst nuclear burning. We note that the small in- 
ferred size of the boundary layer agrees with the theoret- 
ical expectation of most of the boundary layer models at 
sub-Eddington accretion rat es (Kluzniak & Wilson 199l[; 
linogamov fc SunvaevI [19991 : IPopham fc Sunvaev .2001il . 
Thus it is quite possible that the boundary layer emission 
area is indeed small for our observations. 

If the behavior of the hardening factor for the bound- 
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ary layer is similar to that for the burst emission, then 
one might conclude that the boundary layer emission 
area is constant if Lbb Tgg is measured, as the hard- 
ening factor is quite independent of the temperature 
for burst emission in sub-Eddington limit (jMadei et alJ 
120041: lOze]||2006[ ). For example, for a NS with ma,ss 1.4 
Meand radius 10 km. Table 2 from iMadei et all (|2004[ ) 
suggests an approximate relation of Lbb oc Tgg for 
Tbb within 1.1-2.5 keV (most of our observations fall 
into this range). Thus, our results of roughly following 
Lbb oc Tgg from the hard to soft states might imply lit- 
tle change of the real boundary layer emission area for 
our observations. 

We note that there is scatter in the inferred apparent 
emission area of the boundary layer. The fractional vari- 
ation is about 40%, larger than the typical error bars 
(10%). There are also systematic differences in the in- 
ferred boundary layer emission area between BeppoSAX 
and Suzaku, both in the hard and soft states (See Ta- 
bles [3H5]). Whether all this is real or affected by sys- 
tematic problems related to spectral models and/or in- 
strumentation is unclear. If we do a simple power-law 
fit as we did for the MCD component, the BB compo- 
nent seems to follow Lbb cx T^^^^'^ from fit to all data 
combined and Lbb k T^^^^'^ from fit to Suzaku data 
only. Because of the above uncertainties, we did not 
treat these deviations from Lbb c< Tgg as significant. A 
better understanding of our luminosity vs. temperature 
results for the disk and boundary layer of an accreting 
NS may be gained from further theoretical work on each 
component, and from additional observations, e.g., with 
improved statistics and dynamic range for samples of the 
hard state. 

5. CONCLUSION 

The broad-band X-ray spectra of 4U 1705-44 obtained 
with Suzaku and BeppoSAX show two distinct hard and 



soft spectral states. These spectra have significantly bet- 
ter coverage in the soft X-ray energy band compared with 
those from RXTE . We ha ve successfully fit these spectra 
using the model similar to iLin et ahl (|2007| ) . 

The accretion disk in the soft state seems to approxi- 
mately follow a L (X T^'^ track. One cause of the devia- 
tion from L (xT'^ maybe a luminosity-dependent spectral 
hardening factor. However, it is still possible that the 
inner disk radius is really changing in some part of the 
soft state. We found no significant contribution of the 
thermal disk in our hard-state spectra above 1 keV, and 
the disk might be truncated at a large radius and/or has 
a low temperature (< 0.2 keV), or is buried under high 
Comptonization. The boundary layer is roughly constant 
from the hard to soft states, with apparent emission size 
about 1/11 of the whole surface of the neutron star. As- 
suming that the disk is truncated by the ISCO in the soft 
state or the NS surface, we estimated the magnetic field 
of in 4U 1705-44 to be less than about 1.9 x 10^ G. 

Broad relativistic Fe lines are also detected in most 
of the spectra, especially in the soft state. We modeled 
them with the diskline model and found that the strength 
of the Fe line correlates well with the boundary layer 
emission in the soft state, with the Fe line fiux about 
4% of the fiux from the boundary layer (> 6.4 keV). In 
the hard state, our results suggest that the Fe lines are 
due to the strong Comptonization emission. However, 
the Comptonization emission in the hard state seems to 
illuminate the accretion disk and produce the Fe line not 
as efficiently as the boundary layer emission in the soft 
state, probably because the boundary layer is closer to 
the inner accretion disk. 

The authors would like to express their thanks to all 
members of the Suzaku team members, especially Koji 
Mukai, for their support in the schedule of observations 
and preparation of this paper. Support for this research 
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